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ABSTRACT 



Context. Copper is an element whose interesting evolution with metallicity is not fully understood. Observations of copper abundances 
rely on a very limited number of lines, the strongest are the Cu i lines of Mult. 1 at 324.7 nm and 327.3 nm which can be measured 
even at extremely low metallicities. 

Aims. We investigate the quality of these lines as abundance indicators. 

Methods. We measure these lines in two turn-off (TO) stars in the Globular Cluster NGC 6752 and two TO stars in the Globular 
Cluster NGC 6397 and derive abundances with 3D hydrodynamical model atmospheres computed with the C05BOLD code. These 
abundances are compared to the Cu abundances measured in giant stars of the same clusters, using the lines of Mult. 2 at 
510.5nm and 578.2 nm. 

Results. The abundances derived from the lines of Mult. 1 in TO stars differ from the abundances of giants of the same clusters. 
This is true both using C05BOLD models and using traditional 1 D model atmospheres. The LTE 3D corrections for TO stars are 
large, while they are small for giant stars. 

Conclusions. The Cu i resonance lines of Mult. 1 are not reliable abundance indicators. It is likely that departures from LTE should 
be taken into account to properly describe these lines, although it is not clear if these alone can account for the observations. 
An investigation of these departures is indeed encouraged for both dwarfs and giants. Our recommendation to those interested 
in the study of the evolution of copper abundances is to rely on the measurements in giants, based on the lines of Mult. 2. We 
caution, however, that NLTE studies may imply a revision in all the Cu abundances, both in dwarfs and giants. 

Key words. Hydrodynamics - Line: formation - Stars: abundances - Galaxy: globular clusters - NGC 6397, NGC 6752 



1. Introduction 

There is no wide consensus on the nucleosynthetic origin of cop- 
per, and the complex picture drawn by the observations has no 
straightforward interpretation. Multiple channels can contribute 
to the production of this element. According to iBisterzo et alj 
(2004) there are five such channels: explosive nucleosynthesis, 
either in Type II supernovae (SNII) or in Type la supernovae 
(SNIa), slow neutron capture (.s-process), either weak (i.e. tak- 
ing place in massive stars in conditions of hydrostatic equilib- 
rium during He and C burning) or main (i.e. occurring in the 
inter-shell region of low-mass asymptotic giant branch stars) 
and the weak sr-process. The latter occurs in massive stars in 
the C-burning shell when neutron densities reach very high val- 
ues, intermediate b etween typica l s-process neutron densities 
(10 9 - 10 11 cm" 3 : lDespainlll980h and r-pr ocess neutron den- 
sities (10 20 - 10 30 cnT 3 ; iKratz et al.ll2007l) . The contribution 
of the s-process, both we ak and main, to t he sola r system Cu 
abundance is estimated bv iTravaglio et al.l d2004al) to be 27%. 
Explosive nucleosynthesis in SNII can account for 5% to 10% 



Send offprint requests to: Piercarlo.Bonifacio@obspm.fr 

* Based on observations made with the ESO Very Large Telescope at 
Paranal Observatory, Chile (Programmes 71.D-0155, 75.D-0807, 76.B- 
0133) 

** Gliese Fellow 



of the solar system Cu (Bist erzo et al.|[2004h . The contribution 
of SNIa is probably less well known , however, as pointed out by 
iMcWilliam & Smecker-Hand d2005l). the available SNIa yields 
of Cu are rather low dTravaglio et al.ll2004bt iThielemann et al.l 
Il986l) . IBisterzo et all (12004 claim~that the bulk of cosmic Cu 
has indeed been produced by the weak sr-process. 

Observations of copper abundances in Galactic stars show a 
decrease in th e Cu/Fe ratio at low metallicities. This was first 
sugg ested by ICohenl 0980), on the basis of the measurements 
in gian t stars of several Globula r Clusters of different metal- 
licities dCohenl[l978lll979lll980l). It was not until the compre- 
hensive study of ISneden et al.1 (119911) that this trend was clearly 
defined in a robust way, resting on m easurements in a large 
sample of stars. R ecent studies of field (Mi shenina etalj[2002; 
Bihai n et alj |2004|) and Globular Cluster star s ( Shetrone et al.l 
boOlL 120031 ISimmerer et aDl2003L lYong et al.ll2005l) have con- 
firmed this trend (see Fig. 1 of IBisterzo et al.l 20041) . Somewhat 
at odds with these gen eral results are the observations of th e 
Globular Cluster cu Cen (ICunha et al.l2002tlPancino et al.l2002l) . 
Even though this cluster shows a sizea ble spread in metal abun- 
dances (-2.20 < [Fe/H] < -0.70. Ilohnson et alj l200l . the 
Cu/Fe abundance ratio is nearly constant, with no discernible 
trend. 

Observations in Local Group galaxies dShetrone et al.ll200ll 
120031) show that metal-poor populations display low Cu/Fe ra- 
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Fig. 1. Cui 324.7 nm line in the programme stars. The spectra 
are displaced vertically by 0.4 units, with respect to each other, 
for display purposes. 
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Fig. 2. Cui 327.3 nm line in the programme stars. The spectra 
are displaced vertically by 0.4 units, with respect to each other, 
for display purposes. 



tios, similar to what is o bserved in Galactic stars of compara- 
ble metallicity. However, Mc William & Smecker-Hand (120051) 
noted that the metal-rich population of the Sgr dSph displays 
considerably lower Cu/Fe ratios than Galactic stars of compa- 
ra ble metallicity. This re sult is confirmed by the measurements 
of lSbordone etail (120071) . who also include stars of the Globular 
Cluster Terzan 7, associated to the Sgr dSph. 

The majority of the Cu measurements are based on the Cui 
lines of Mult. ^J, sometimes one line of Mult. 7 i s used . The 
exceptions are the measurements of Bih ain et alj |2004) and 
Cohen etaD, d2008|). w ho use the resonance lines of Mult. 1 and 
Prochaska et al. (2000) who, to our knowledge, are the only ones 
who have made use of the strongest line of Mult. 6 in the near in- 
frared. While for stars of metallicity above -1.0 one may have a 
choice of several lines to use, when going to metal-poor stars, 
for instance below -1.5, the only usable Cu abundance indica- 
tors are the two strongest lines of Cu i Mult. 2 in giant stars and 
the resonance lines of Mult. 1 in both dwarfs and giants. The 
advantage of the lines of Mult. 1 is that they are very strong, 
at high metallicity they are strongly saturated and therefore not 
ideal for abundance work, but, they remain measurable down to 



an extremely low metallicity. Bihain et al. (2004) have been able 
to measure the 327.3 nm line in the extremely metal-poor dwarf 
G64-12 ([Fe/H]~ -3). Observationally the main disadvantage 
of Mult. 1 is that it lies in the UV, fairly near to the atmospheric 
cut-off. A very efficient UV spectrograph, like UVES and a large 
telescope, like VLT, may circumvent this problem. There are 
many spectra suitable for the measurement of the Cu i lines of 
Mult. 1 in stars of different metallicities in the ESO archival 
The main purpose of our investigation is to assess the quality of 
the Cui lines of Mult. 1 as abundance indicators. Our strategy is 
to compare for the first time Cu abundances in main sequence 
and giants of the same cluster, because Cu is not expected to be 
easily destroyed or created. This test will indicate the reliability 
of our modelling. Globular Clusters NGC 6397 and NGC 6752 
span an interesting range in metallicity -2.0 to -1.5, which is 
relevant for a large fraction of the observations in field stars. 



we refer to the multiplet designation of Moore ( 1945) 
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Table 1. Atmospheric parameters of the programme stars. 



Star 


T 

1 eff 




TFe/Hl 
[rc/nj 


P 




K 


[cgs] 


dex 


kms~' 


CI* NGC 6752 GVS 4428 


6226 


4.28 


-1.52 


0.70 


CI* NGC 6752 GVS 200613 


6226 


4.28 


-1.56 


0.70 


CI* NGC 6397 ALA 1406 


6345 


4.10 


-2.05 


1.32 


CI* NGC 6397 ALA 228 


6274 


4.10 


-2.05 


1.32 


CI* NGC 6397 ALA 2111 


6207 


4.10 


-2.01 


1.32 


HD 218502 


6296 


4.13 


-1.85 


1.00 



2. Observations and equivalent width 
measurements 

The sp ectra analysed here are th e same as in iPasquini et al.l 
(2004) and Pasquini et al. (2007) and described in the above 
papers. They were obtained with the UVES spectrograph 
ilDekkeret aLll2000h at the ESO VLT-Kueyen 8.2m telescope. 
We here use the blue arm spectra, which are centred at 346 nm. 
Both clusters were observed with a 1" slit and a 2 x 2 on-chip 
binning, which yields a resolution of about 40 000. The reduced 
spectra were downloaded from the ESO archive, thanks to th e 
improved strategies for optimal extraction dBallester et al.l2 006). 
the S/N ratios are greatly improved compared to what was pre- 
viously available. The equivalent widths (EWs) of the two Cu i 
lines of Mult. 1 were measured with the IRAF task splot and 
are provided in the on-line Table [6] 

In addition to the four cluster stars we analyse the field star 
HD 218502 as a reference. Its atmospheric parameters are close 
to those of the cluster stars. For this star we work with the 
data used by Pasqu ini et al.1 (120041) as well as with the data ob- 
served in 2005, in th e course of ESO programme 76.B-0133 (see 
Smili anic et al.l2008h . For this star we used six spectra: two with 
l'.'O and 2x2 binning, two with l'.'O and lxl binning, two with 
1'.'2 and lxl binning. Each pair of spectra was coadded, the 
equivalent widths were measured on the coadded spectrum and 
then the three equivalent widths were averaged. The spectra of 
all the five stars analysed here are shown in Fig.[TJand[2] 

One of the goals of the present analysis is to compare the Cu 
abundances in the TO stars with those measured in giant stars of 
the sa me cluster. For NG C 6752 we can rely on the recent anal- 
ysis bv lYong et al.l (120051) . who analysed 38 giants in this clus- 
ter, making use of high-reso lution high-S/N ra tio UVES spec- 
tra. The atomic data used by I Yong etail (12001 are the same as 
those here used. The analysis is based on ID ATLAS models 
an d LTE spect r um sy nthesis. That the ATLAS models employed 
by lYong et aTl (120051) use the approximate overshooting option 
in ATLAS, while those we use do not , brings about only minor 
differences. Thus the measurements of Yon g et al.1 (120051) are di- 
rectly comparable to our own. The measurements of lYong et al.l 
(120051) are, however, based only on the strongest line of Mult. 2. 
Because we aim to compare the abundance derived from the 
lines of Mult. 1 and Mult. 2 we retrieved UVES reduced spectra 
of one of the stars o f lYong et al.l J2005) from the ESO archive: 
CI* NGC 6752 YGN 30. We used three spectra of 1800 s ob- 
tained with the dichroic # 1 , the blue arm spectrum was centred 
at 346 nm and the red arm spectrum at 580 nm. The slit was set at 
l'.'O in the blue arm and (X'7 in the red arm; the CCD binning was 
lxl for both arms. The corresponding resolution is ~ 45 000 
for the blue arm and ~ 60 000 in the red arm. 

For the cluster NGC 6397, though, we were unable to find 
any recent analysis that included the measurement of Cu. In fact 
the only measurement of Cu in this cluster which we could find 



is due to iGrattonl (Il982l) . In order to make the iGrattonl (1 19821) 
measurements directly comparable to our own we used the pub- 
lished EWs of the Cu i lines of Mult. 2 and derive the abundances 
with our models, spectrum synthesis codes, and atomic data. 

3. Cu abundances 

3.1. Atomic data 

To determine the Cu abundances for the TO stars we used the 
Cui resonance lines of Mult. 1 at 324.7 nm and 327.3 nm. 
The log gf values were taken from Bielski (1 19751) and the hy- 
perfine stru cture and iso topic shifts for the 63 Cu and 65 Cu iso- 
topes from iKurucz I dl999l) . We used the same sources for the 
two lines of Mult. 2 at 510.5 nm and 587.2 nm that we used for 
the giant stars. The line list used for the computations is given 
on-line in Table [5] The line at 327.3 nm is free from blends 
in metal-poor TO stars and the continuum is usually easily de- 
termined. The stronger 324.7 nm line, though, lies in a more 
complex spectral region. The only truly blending feature is a 
weak OH line (324.761 5 nm), but, the line is on the red wing 
of a complex blend, mainly of iron lines, of which several have 
poor log gf values. The continuum is more difficult to deter- 
mine for this line, given the larger line crowding in this region. 
We experimented with different choices for the Van der Waals 
broadening of the lines, the ABO theory dAnstee &Q' Mara 
ll995blBarklem & 0'Marall997l:|Barklem et alJl998allbl) and the 
WI DTH approxim ation dKuruc d 1 993al 120051: ICastellill2005L see 
also lRvanl ll998). For the transitions under consideration the 
WIDTH approximation and the ABO theory yield almost iden- 
tical values. 

3.2. Atmospheric parameters 

The adopted atmospheric parameters for our programme stars 
are given in Table p] and were taken from Pasqui ni et al.l (2004) 
and IPasquini et alfihoOTl) . For the giant star CI* NGC 6752 
YGN 30 we adopted the atmospheric parameters of Yon g et all 
d2005l) . For the two gi ants in NGC 63 97 we adopted the atmo- 
spheric parameters of Grattqn] (Il982l) . For the reader's conve- 
nience the atmospheric parameters of the giant stars are provided 
here on-line in Table [7] 

3.3. Model atmospheres and spectrum synthesis. 

For each star we computed a ID model atmosp here using ver- 
sion 9 of the ATLAS code dKuruczlll993al 120051) under Linux 
dSbordone et al.ll2004l:ISbordonel2005l). We used the opacity dis- 
tribution functions described by ICastelli & Kuruczl (120031) and 
microturbulent velocity 1 kms -1 , the mixing-length parame- 
ter, aMLT, was set to 1.25, and the overshooting switched off. 
This model atmospher e was used as input to the SYNTHE code 
dKuruczlll993bl 120051) . with different Cu abundances, to com- 
pute a curve-of-growth for each line. The Cu abundances were 
derived by interpolating in these curves of growth. The corre- 
sponding abundances are given in the second column of Table 
[2] the cr is the variance of the abundances of the two lines. The 
abundances for the individual lines can be found on-line in Table 
i 

The use of three dimensional hydrodynamical simula- 
tions to describe stellar atmospheres (hereafter 3D models) 
has led to the important notion that the outer layers present 
steeper temperature gradie nts than predicted by traditional ID 
static model atmospheres ( Asplu nd et alJ[l999l : lAsplundl f2005: 
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Table 2. Copper abundances for the programme stars. 
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Fig. 3. Curves of growth (COG) for the Cui 327.3 nm transi- 
tion. The dotted line to the left is the COG for the 3D model 
d3t63g40m20n01, the three solid lines to the right are those for 
the corresponding IDlhd model for three values of the micro- 
turbulent velocity: 0.5, 1.0 and 1.5 kms -1 , bottom to top. The 
dashed line is the 3D COG shifted arbitrarily by +0.58 dex along 
the x-axis. This highlights that the shape of the 3D COG differs 
from that of the corresponding ID COGs and therefore the 3D 
correction depends on the EW of the transition. 

iGonzalez Hernandez et a and that this effect is consid- 

erably more pronounced for metal-poor stars. In addition to the 
different mean temperature profile, the 3D models differ from 
traditional ID models because they account for the horizontal 
temperature fluctuations. Both effects may or may not be impor- 
tant, depending on the line formation properties of the transition 
under consideration. In order to investigate these effects for the 
Cui lines w e used several 3D models computed with the cod e 
C0 5 BOLD dFrevtag et al.ll2002l 120031: fWedemever et al.ll2004h . 
The characteristics of the 3D models employed in this study are 
given in Table [3] The line formation computations for the 3D 
models were performed with the Linfor3D code0. For each 3D 
model we used also two reference ID models: the (3D) and the 
IDlhd, which we define below. 

The (3D) models are computed on-the-fly by Linfor3D by av- 
eraging the 3D model over surfaces of equal Rosseland optical 
depth and time. The (3D) model has, by construction, the mean 
temperature structure of the CO 5 BOLD model, therefore the dif- 
ference in abundance A(3D)-A((3D)), allows us to single out the 
effect s caused by temperature fluctuations (see Caffau & Ludwig 
120071) . 

The IDlhd model is a ID, plane parallel, LTE, static, model 
atmosphere and employs the same micro-physics and opacity 
as the CO 5 BOLD models; it is computed with the LHD code. 
These models are our models of choice to define the "3D cor- 
rection" as A(3D) - A(1DlhdX where A is the abundance of any 
gi ven element. More detail s on the LHD models m ay be found 
in lCaffau & Ludwigl d2007h and lCaffau et alj (12010b . 

In any given Linfor3D run we made computations also for the 
(3D) model and for a IDlhd model, with the same T e jf, logg, 
and metallicity as the 3D model. 

The computation of a 3D model is still very time consum- 
ing, even on modern computers (several months), it would be 
impractical to compute a specific 3D model for any set of our 

3 http://www.aip.de/~mst/Linfor3D/linfor_3D_manual.pdf 



Star 




A(Cu) 


(T 




(X 






ID 


3D 




CI* NGC 6752 GVS 4428 


3.23 


0.08 


2.56 


0.16 


CI* NGC 6752 GVS 200613 


3.01 


0.05 


2.23 


0.07 


CI* NGC 6397 ALA 1406 


1.33 


0.03 


0.74 


0.05 


CI* NGC 6397 ALA 228 


1.30 


0.03 


0.73 


0.05 


CI* NGC 6397 ALA 21 11 


1.19 


0.02 


0.60 


U.U2 


HD 218502 




1.52 


0.09 


0.95 


0.04 


Table 4. Mean copper abundances for the two clusters 




Star 


A(Cu) 


a 


A(Cu) 


a 






ID 




3D 






NGC 6752 dwarfs 


3.04 


0.07 


2.28 


0.12 




NGC 6752 giants 


2.03 


0.05 


1.98 


0.05 




NGC 6397 dwarfs 


1.25 


0.05 


0.63 


0.04 




NGC 6397 giants 


1.40 


0.17 


1.30 


0.17 





atmospheric parameters. Our strategy is therefore the following: 
we perform an abundance analysis with ATLAS model atmo- 
spheres computed for the desired set of atmospheric parameters; 
we use a grid of 3D models with atmospheric parameters that 
bracket the desired ones and compute the relevant 3D correc- 
tions by linear or bi-linear interpolation in the grid, as appropri- 
ate; the 3D abundance is obtained by applying the interpolated 
3D correction to the ID abundance. 

The line formation computations were performed using 
SYNTHE for the ATLAS models and Linfor3D for all other mod- 
els. As a consistency check we used Linfor3D with an ATLAS 
model as input and verified that the line profiles and EWs are 
consistent with the results derived from SYNTHE+ATLAS. The 
difference between the two line formation codes amounts to a 
few hundredths of dex in terms of abundance, a quantity that is 
irrelevant with respect to the size of the 3D corrections under 
consideration. 

For three out of the six TO stars under study, the Cu 1 lines 
are strong (EW> 4.0pm), therefore they are surely in the satura- 
tion regime. Their 3D correction depends on the adopted micro- 
turbulence in the adopted reference ID atmosphere. To take this 
into account, different IDlhd curves of growth were computed 
with microturbulent velocities of 0.5, 1.0 and 1.5 kms -1 , to al- 
low us interpolation to any desired value of For weaker lines 
the microturbulent velocity does not play a fundamental rule, so 
that the 3D correction is mostly insensitive on the choice of this 
parameter. We find that the 3D curve of growth is not a simple 
translation of a ID curve of growth, but has a distinct shape. An 
example to illustrate the effect is shown in Fig. [3] An immediate 
consequence is that for all lines we considered, the 3D correc- 
tion depends on the EW, even for the weaker lines. We note that 
the dependence of the 3D correction on the microturbulence im- 
plies that the abundance obtained by applying the correction to 
an abundance derived from a ID model depends on the adopted 
microturbulence. One of the reasons to prefer the use of 3D mod- 
els is to avoid the use of this parameter. The correct way to treat 
this is, not to use the 3D correction, but to derive the abundance 
by an interpolation in a set of 3D curves of gr owth, or to use 
suitab le fitting functions, as done for lithium by Sbor done et al.l 
(l2010h . But this requires the use of a larger set of 3D models, 
which bracket the effective temperatures of the studied stars. For 
the purpose of the present exploratory investigation we believe 
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Table 3. CO 5 BOLD models employed in the study. 



Model 




logg 


[M/H] 


N t 


time 


tc 


Resolution 


Box Size 




K 








s 


S 




Mm 3 


d3t50g25mmlOn01 


4990 


2.5 


-1.0 


20 


475990 


1411.9 


160 x 160 x 200 


573.2 x 573.2 x 245.4 


d3t50g25mm20n01 


5020 


2.5 


-2.0 


20 


403990 


1388.3 


160 x 160 x 200 


584.0 x 584.0 x 245.4 


d3t63g40mml0n01 


6260 


4.0 


-1.0 


20 


43800 


12.2 


140 x 140 x 150 


26.0 x 26.0 x 12.8 


d3t63g40mm20n01 


6280 


4.0 


-2.0 


16 


27600 


49.0 


140 x 140 x 150 


26.1 x 26.1 x 12.8 


d3t63g45mml0n01 


6240 


4.5 


-1.0 


20 


24960 


16.0 


140 x 140 x 150 


7.0x7.0x4.0 


d3t63g45mm20n01 


6320 


4.5 


-2.0 


19 


9120 


15.9 


140 x 140 x 150 


7.0x7.0x4.0 



the approach of using 3D corrections is adequate, especially be- 
cause the conclusion of the study is that 3D-LTE abundances 
are unreliable. Finally we point out that it may be that current 
3D models do not correctly capture turbulence at small scales 
(Steff en et al] l2009). If this were indeed the case all abundances 
derived from saturated lines are doubtful, whether they are de- 
rived applying a 3D correction or directly derived from the 3D 
curves of growth. 

We therefore computed the 3D correction for each of the 
measured EWs for each of the relevant 3D models. In general the 
3D correction will also depend on the treatment of convection in 
the ID reference model, hence on the adopted ffMLT- The lines 
under consideration do not form in the deepest layers, which are 
the most affected by the choice of Q'mlt? thus they are insensitive 
to it. All IDlhd models employed have aMLT= 1-0. The com- 
puted A(3D) - A( IDlhd) corrections, as well as the A(3D) - 
A((3D)) corrections are given on-line in Table [8] The 3D abun- 
dances provided in Table [2] are obtained by applying to the ID 
abundances the 3D corrections in Table [6] which were obtained 
by interpolating the corrections in Table [8] 

The six stars under study have very similar effective temper- 
atures. All are within roughly 100 K of the effective temperatures 
of the 3D models listed in Table[3](T eff ~ 6300 K). Therefore it 
is not necessary to include more 3D models and perform an in- 
terpolation in T e ff . The metallicities and gravities of the models 
in Table [3] bracket the metallicities and gravities in Table [TJ We 
used a bi-linear interpolation in metallicity and gravity for the 
two stars of NGC 6752 and for HD 218502. The three stars of 
NGC 6397, though have a metallicity of almost -2.0, therefore a 
linear interpolation in surface gravity was sufficient. 

The computation of 3D models for typical giant stars is much 
more time consuming than for F-type and cooler dwarfs. The ra- 
diative relaxation time in the surface layers of warm giants be- 
comes significantly shorter than the dynamical time scale, which 
makes it computationally expensive to properly capture the time 
evolution of the system. At present we have only two fully re- 
laxed models of metal-poor giants. The parameters are given in 
Tableland the metallicities are -1 .0 and -2.0. The surface grav- 
ity is larger than that of the majority of the giant star s analysed 
in eith er cluster. Only the faintest stars analysed by lYong et al] 
(2005) have parameters T e ff and logg close to those of our gi- 
ant models. In spite of this we believe that we can use the 3D 
corrections derived from this model, as a first order approxima- 
tion as representative of the corrections in giant stars of both 
clusters. This is possible because the 3D corrections are rather 
small, especially when compared to those of dwarf stars. This 
is partly because giant 3D models do not show a pronounced 
over-cooling, compared to ID models, as dwarfs display. This 
conclusion is also based on the examination of several snapshots 



E 




Fig. 4. Contribution functions of the EW at disc-centre, de- 
fined in a way that their integral over log ta gives the EW 
(Magain 1986), for the Cui 324.7 nm line and the model 
d3t63g40mm20n01 for two different values of Cu abundance. 
In the upper panel A(Cu)=0.2, in the lower panel A(Cu)=1.7. 
The solid lines refer to the 3D model, the dashed lines to the 
corresponding IDlhd model. 
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Fig. 5. Contribution functions for the Cu i 5 10.5 nm line and the 
model d3t50g25mm20n01. The solid line refers to the 3D model, 
the dashed line to the corresponding IDlhd model. 



of not fully relaxed giants of different atmospheric parameters, 
which we are in the process of computing. The 3D correction for 
the giant models is -0.1 dex for both the examined Cui lines of 
Mult. 2 for the model of metallicity -2.0 and 0.0 dex (actually 
+0.01) for the model of metallicity -1.0. We apply a correction 
of -0. 1 to the abundances of giant stars in NGC 6397 and -0.05 
to those of NGC 6752. 
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4. Results 

4.1. HD 218502 

The analysis of the reference star HD 218502 shows that the 
abundances derived from the two Cui resonance lines (Table [6} 
agree well, both with ID and 3D models. This gives us con- 
fidence in the reliability of the atomic data used. It also sug- 
gests that with good quality data, the EWs of both lines can be 
satisfactorily measured in spite of the complexity of the spec- 
tral region, especially for the 324.7 nm line. Our ID abundance 
(A (Cu)= 1.52 + 0.09) agrees, within errors with what reported 
bv lBihain et all (12004 A(Cu)= 1 .70 + 0. 1 7). We note a small dif- 
ference in the effective temperature adopted in the two analyses 
(about 100 K) and a difference in the data: our UVES spectra are 
of considerably hig her quality than the CASPEC spectra used by 
Bihai n et alj ([2004), who measured only the 327.3 nm line. This 
agreement is expected because the atomic data are the same in 
the two analyses and the ID model atmospheres used are simi- 
lar, the main difference being the overshooting. This check sug- 
gests that we may reasonably assu me that our Cu ab undances 
should be consistent with those of Bihain et al . (2004), which 
are based on the UV resonance lines. An inspection of Fig. 6 of 
Bihai n et al.l (|2004) sugg ests that these measure ments substan- 
tially agree with those of Mishenina et al ] (120021) . which are es- 
sentially based on the measurements of the lines of Mult. 2. Yet 
one should take into account the large error bars that derive from 
the relatively poor S/N ratio that is achievable in the UV range. 

The 3D correction for the lines of Mult.l is large and the 
3D abundance in this star is well below all the measurements 
in giant stars of similar metal licity. An applicati on of 3D correc- 
tions to all measurements of Bihai n et all d2004l) would probably 
break the agreement with the measurements of [Mishenina et al] 
(2002). We have inspected the available red spectra of HD 
218502, to see if any of the lines of Mult. 2 could be detected, 
but this was not the case. 

4.2. TO stars in NGC 6752 and NGC 6397 

For the cluster stars there is also a good consistency between the 
abundances derived from the two resonance lines for any given 
line, in spite of the much lower S/N ratios in the cluster stars. 
This suggests that there is no major inconsistency in the EW 
measurements. 

The weighted mearQ of the abundances of the clusters, re- 
ported in TableHJ displays an error in the mean that is reasonably 
small. We believe that the mean abundances for the two clusters 
obtained from the dwarf stars are indeed representative of the 
Cu i abundance derived from the lines of Mult. 1 . Considering 
that each spectrum of a TO star in these clusters amounts to 
about 10 hours of integration with UVES, it is unlikely that in 
the near future better data or data for a larger number of stars 
will be available, although it would clearly be desirable. 

4.3. Star CI* NGC 6752 YGN 30 

In this star, like in the other giant stars observed by Yon g et aTl 
(2005), both the lines of Mult. 1 and of Mult. 2 are measurable, 
which provides for a consistency check. We did not measure the 

4 The error on the weighted mean has been taken to be the largest 

between and ^r) x ^rj) 2 ^p^- where ** 

the data points cr, are the associated e rrors, (x) is th e mean value and n 
is the number of data point (see e.g., [Agekan 1972, pages 144—150). 



line at 324.7 nm since the region is extremely crowded in these 
cool giants, but the line at 327.3 nm is clean and unblended. For 
Mult. 2 we only measured the 510.5 nm line, since the other line 
is not present in the spectrum, because it falls in the gap between 
the two CCDs. The two lines (324.7 nm and 5 10.5 nm ) provide 
inconsistent results, the line of Mult. 1 provides an abundance 
that is 0.54 dex higher than that of the line of Mult. 2 (see on- 
line Table The abundance we derive from the line of M ult. 2 
substantially agrees with the measure of lYong et al.l d2005l) . our 
abundance is 0.13 dex smaller. Of this difference 0.05 dex are 
because w e use ATLAS models without overshooting, while 
Yon g et al.l (120051) use models with overshooting, the remaining 
difference should be attributed to a difference in the measured 
EW and possib ly to the different spectrum synthesis codes used 
(MOOG bv lYongetal]|2005l and SYNTHE by us). The 3D cor- 
rections for both lines of Mult. 1 and Mult. 2 are small and agree 
to within a few hundredths of dex. The abundance from the two 
lines cannot be brought into agreement by using 3D models. 

4.4. Giant stars In NGC 6397 

The error in the mean of the two gia nts is . 18 dex , which is 
essentially identical to the estimate of Gratton (1982) of 0.2 dex 
on the copper abundance. Even making use of ID models the 
abundance of the giant stars is considerably higher than in the 
dwarf stars. 



5. Effects of atmospheric parameters 

We intend to quantify the effect of changing atmospheric param- 
eters on the derived abundances. The abundances derived from 
the Cu i resonance lines are fairly sensitive, for a neutral species, 
to the adopted surface gravity. For dwarf stars a change of +0.25 
dex in log g induces a change of +0.1 dex in abundance. For the 
giant stars they are only slightly less sensitive, +0.06 dex. On the 
other hand for giant stars the dependence on gravity of the abun- 
dances derived from the lines of Mult. 2 is very weak +0.01 dex. 
An inspection of on-line Table [8] allows us to estimate the effect 
of surface gravity on the 3D corrections for the lines of Mult. 1 . 
By increasing the gravity by 0.5 dex, the 3D correction increases 
by 0.1 to 0.2 dex, depending on how saturated the line is. Since 
the 3D correction is negative, this means that it decreases in ab- 
solute value. The opposite trends with surface gravity on the ID 
abundance and 3D correction imply that the two effects tend to 
cancel and the overall sensitivity of the 3D abundance on surface 
gravity is small. 

The dependence of abundances on effective temperatures for 
the lines of Mult. 1 is similar for dwarfs and giants, and is about 
+0.2 dex for a change of +100 K in effective temperature. To 
evaluate the dependence of the 3D corrections on the effective 
temperature we used four models, extracted from the CIFIST 

I . 1 1 

grid (Ludwig et al. 2009). All four models have effective tem- 
perature around 5900 K, their metallicities are -1.0 and -2.0 and 
their log g 4.0 and 4.5. The result is that for a decrease of 300 K 
the 3D correction increases by 0.2 dex. Again the variation of the 
3D correction goes in the opposite direction with respect to the 
variation in the ID abundance. Combining the results we con- 
clude that a decrease of 300 K in effective temperature results in 
a decrease by 0.4 dex in copper abundance (3D). 

For the giants stars we also estimated the variation of the 
abundances derived from the lines of Mult. 2, which amounts 
to about +0.1 dex for a variation of ±100 K in effective tempera- 
ture. 
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T etf = 6300 log g = 4.0 [M/H] = -2 
12 bin model 



-6 -4 -2 



Fig. 6. Comparison between the temperature structures of two 
models with T e ff — 6300, log g= 4.0 and metallicity -2.0, com- 
puted with 6 and 12 opacity difference. The differ ence is not as 
large as it is at lower metallicities (see Fig. 1 of iBehara et al.l 
l2009h . 

6. Hydrodynamical models and spectrum 
synthesis. 

Clearly the large 3D corrections derived for the Cu i resonance 
lines are driven by the fact that in the outer layers the hydrody- 
namical models are considerably cooler than the corresponding 
ID models, the so-called "overcooling". This shifts the ionisa- 
tion equilibrium towards neutral copper, and as a consequence 
the line contribution function shows a strong peak in these lay- 
ers, contrary to the ID model (se Fig.|4j». From a physical point 
of view this is expected, simply because the hydrodynamical 
model is observed to transport flux through convection even 
in layers where the corresponding ID model is formally sta- 
ble against convection (overshooting). One should however ask 
to which extent the compu t ed ove rcooling depends on the as- 
sumptions made. Bonifacio (2010) pointed out the difference in 
the overcooling for met al-poor gian t s com p uted with CO 5 BOLD 
and that computed by ICollet et all (120071) . IBehara et all (120091) 
pointed out that for extremely metal-poor dwarfs the overcool- 
ing is considerably less in CO 5 BOLD models computed with 12 
opacity bins than in models computed with 6 opacity bins, like 
the ones used here. 

For the time being we have relatively few CO 5 BOLD models 
computed with 12 opacity bins, among which we have one with 
parameters identical to one used in the present investigation: 
T e ff= 6300, logg= 4.0 and metallicity = -2.0. In Fig.|6]we show 
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Fig. 7. Comparison between the curves of growth for the 
327.3 nm line computed for the two models shown in Fig. [6] 
Solid symbols correspond to the 6 bin model and open symbols 
to the 1 2 bin model. 



the mean temperature structures of the two models. Obviously 
the difference is smaller than w hat is displayed by the models 
1 dex more metal-poor shown bv lBehara et alj d2009f) . The qual- 
itative conclusion is confirmed by comparison of the curves of 
growth. In Fig.|7]the curve of growth for the 327.3 nm line used 
in the present investigation is compared to the one computed 
from the corresponding 12 bin model. As can be appreciated 
from the plot, for a given equivalent width the 12 bin model will 
yield a Cu abundance that is higher by approximately 0. 1 dex, 
thus correspondingly decreasing the 3D correction. 

Another matter of concern is that the current version of 
Linfor3D treats scattering as true absorption. That is to say, al- 
though scattering processes such as Rayleigh scattering off hy- 
drogen atoms are taken into account as opacity sources, in the 
solution of the transfer equation the source function is set equal 
to the local Planck function, without a term depending on the 
mean radiation intensity (5 y = B Y ). To which extent can this ap- 
proximation affect our computations, especially in the near ultra- 
violet, where scattering processes are a non-negligible source of 
opacity ? We assessed this by using ID models and ID spectrum 
synthesis. We used a slightly modified version of the SPECTRV 
code in the SYNTHE suite so that scattering is treated as true 
absorption whe n the card SCA TTERING OFF is set in the input 
model (see e.g. lCastellilll988h . We computed line profiles both 
with SCATTERING ON and SCATTERING OFF for the model 
at T eff = 6296 K logg= 4.0 and metallicity -2.0 (rlevant to HD 
218502) and for the model T eff = 4943 log g= 2.42 and metallic- 
ity -1.5 (relevant to CI* NGC 6752 YGN 30). It turns out that 
in both cases the difference is irrelevant to our analysis (0.6% in 
the continuum and 0.03% in the residual intensity for the dwarf 
model and 14% in the continuum and 2% in the residual inten- 
sity). The effect of treating scattering as true absorption is very 
similar on the continuum and in the lines, thus implying a small 
effect on the residual intensity an d equivalent width. 

Recently Haye k et al.l d2010b have introduced a proper treat- 
ment of scattering in their magneto-hydrodynamical simulation 
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code, BIFROST. For the Sun and solar-type stars they do not find 
a significant impact of continuum scattering on the temperature 
structure. We thus believe that although the impact of a proper 
treatment of scattering needs to be investigated in metal-poor 
dwarfs, it seems unlikely that the results presented here will be 
seriously challenged by this. 

7. Discussion 

If we consider the abundances in Table|4]at face value we are led 
to the inescapable conclusion that the Cu abundances in dwarfs 
and giants do not agree. Even though the abundances are almost 
compatible, within errors, at least in the ID case, this system- 
atic difference should not be overlooked. The behaviour is dif- 
ferent between the two clusters: in NGC 6752 the dwarfs pro- 
vide a higher abundance than the giants, while in the case of 
NGC 6397, they provide a lower value than the giants; this both 
using ID and 3D models. For NGC 6752 the difference in ID 
abundances is of one order of magnitude (1 dex), but this is re- 
duced to only 0.3 dex if we look at the 3D abundances, which 
are almost compatible with the errors. The situation is reversed 
in NGC 6397, in ID the abundances of dwarfs is only 0.15 dex 
lower than that in giants, while in 3D, the diferece is 0.9 dex. 
This behaviour may be understood in terms of the different line 
formation properties in dwarfs and giants and how they change 
with different Cu abundance. In Fig.|4]we show, as an example, 
the contribution functions of the EW at disc-centre for one of 
our models for a dwarf star, for two different Cu abundances. In 
the top panel A(Cu)=0.2 and the line is on the linear part of the 
curve of growth, in the bottom panel, A(Cu)=1.7 and the line 
is saturated. In both cases the 3D contribution function is very 
different from the ID one and is peaked in the outer layers of 
the atmosphere. The formation of the lines of Mult. 2 in the at- 
mospheres of giants is instead much less affected by 3D effects, 
as depicted in Fig. [5] The contribution functions of the lines of 
Mult. 1 in the giant model are morphologically similar to what is 
shown in Fig. [5] confirming the weak overcooling present in this 
model. While the above arguments explain the behaviour of the 
3D corrections, they have no bearing in the abundance difference 
that we find between dwarfs and giants. 

In these situations one has always to consider two possible 
alternatives: i) the difference is true and has an astrophysical ori- 
gin; ii) the difference arises from shortcomings in the analysis. In 
our opinion hypothesis i) must be discarded. It seems extremely 
far fetched to devise a physical mechanism by which Cu should 
be overabundant in dwarf stars, like in NGC 6752, while the re- 
verse is true in NGC 6397. Where giants have a higher Cu abun- 
dance, one could imagine to explain such a scenario either by 
invoking diffusion in TO stars, or Cu production in giant stars, 
or perhaps even a combination of both. However, an examina- 
tion of Table [4] shows that the differences that need explanation 
are far too large to be created by diffusion, and the Cu produc- 
tion would also have to be highly efficien t. In addition note that 
the Cu abundances in the large sample of lYong e t al. (2005|) are 
extremely uniform, which speaks against Cu production. 

Thus we are left with the conclusion that the abundance 
determinations in either dwarfs or giants, or both, are wrong. 
Let us start by examining the 3D abundances in dwarf stars. 
Contribution functions like those shown in Fig.|4]must give rise 
to concerns about the LTE approximation used in our computa- 
tions. The outer and less dense layers of the atmosphere, which 
contribute mostly to the line EW, are those in which the photon 
mean free path is longest and deviations from LTE may be ex- 
pected. The situation is even more extreme in a 3D atmosphere 



where photons from a hot up-draft may transfer horizontally and 
overionise a neighbouring cool down-draft. A morphologically 
similar situ ation is indeed observed for the Li i dou blet in metal- 
poor stars ( Asp lund et alll2003t ICavrel et al]l2007l) . In LTE the 
contribution function displays a double peak, with a substantial 
contribution from outer cool layers, in NLTE this peak is entirely 
suppressed by overionisation. The nearly exact cancellation be- 
tween 3D and NLTE correction that takes place for the Li i dou- 
blet, and results in 3D-NLTE abundances in very close agree- 
ment with 1D-LTE abundances, should not be taken as a general 
rule. Nevertheless we believe that the Cu i lines of Mult. 1 cannot 
be described by 3D-LTE computations, but NLTE effects should 
be properly accounted for. 

This leads to the question of how reliable the LTE approx- 
imation is for the ID computations. That even in ID the abun- 
dances in dwarfs and giants differ by a factor of the order of 2, 
surely prompts to see if for NLTE computat ions the two sets of 
abundances may be brought into agreement. Bih ain et alj ([2004) 
tried to fit the Cu i lines of Mult. 1 in the solar spectrum, but were 
unable to reproduce the core of the lines. This was attributed 
mainly to the presence of a chromosphere that influences the 
cores of strong lines. Deviations from LTE, however, could be 
another, possibly concomitant, cause for the failure to reproduce 
the line core in LTE. For Pop II dwarf stars chromospheric ef- 
fects should not be strong, in view of their old age, even if chro- 
mospheres were present. 

Let us finally consider the Cu abundances in giant stars. Our 
limited computations suggest that they should not suffer large 
3D effects. The test we conducted for the giant star in NGC 
6752 suggests that the LTE synthesis does not allow us to re- 
produce satisfactorily the lines of Mult. 1 and of Mult. 2 with the 
same Cu abundance. Indeed the discrepancy is quite large (0.5 
dex); significant deviations from LTE for either or both sets of 
lines could be responsible for this. Which should be further in- 
vestigated in order to produce reliable Cu abundances. The evo- 
lution of copper with metallicity, essentially based on the mea- 
surements of Mult. 2 in giant stars, shows a rather sharp drop 
in [Cu/Fe] around [Fe/H]=-1.5. This means that it takes place 
around A(Cu)=2.0. In the curve of growth for the 510.5 nm line 
of Mult. 2 for our giant model, this is roughly the abundance for 
which the line begins to enter in the saturation regime. If the 
Cu i lines of Mult. 2 suffer deviations from LTE it is likely that 
these depend on the line strength and may show a rather sharp 
change just when the line enters a saturation regime. This be- 
haviour i s observed, for instance , for the sodium D lines (see 
Fig. 6 of Andrievs kv et alJ l2007). These considerations render 
NLTE computations for Cu very desirable. 

Unfortunately, to our knowledge, up to now no such compu- 
tations have been published, nor does a Cu model atom exist. 

Among the possible causes for the discrepancy in abun- 
dances between giants and dwarfs one may also consider errors 
in the atmospheric parameters. One may conclude that this can- 
not be the case by noticing the discrepancy between the Cu abun- 
dance derived from Mult. 1 and Mult. 2 in star CI* NGC 6752 
YGN 30, for which lines of both multiplets are measured. Given 
that the response of both multiplets to a change in effective tem- 
perature is similar, the conflicting results cannot be resolved by 
changing the effective temperature. Of course the discrepancy 
between the two multiplets can be resolved in ID by invoking 
a higher microtrubulence, but it would be necessary to raise it 
by 1 kms~'. This increase would then cause a strong trend be- 
tween iron abundances and equivalent widths. Furthermore this 
would not allow us to solve the discrepancy in the 3D analysis. 
Indeed while the 3D correction for Mult. 2 would not change sig- 
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nificantly with this increase in microturbulence, those of Mult. 1 
would increase by about 0.4 dex, thus breaking the agreement 
between the two multiple ts forced in ID. Formally one can cer- 
tainly find a value of the microturbulence that forces the ID 
abundance plus 3D correction of the two mutliplets to be equal, 
while leaving a discrepancy in the ID abundances. This however 
would again cause an abundance spread among lines of other el- 
ements (e.g. iron) and can hardly be invoked as a solution of the 
problem. Although at the moment we do not have enough 3D 
models for giant stars to perform a full 3D analysis, as done by 
ISbordone et al.l d2010l) for lithium, we believe that our results in- 
dicate that this analysis will provide discrepant abundances from 
the two multiplets. 

Let us further consider if changes in the atmospheric param- 
eters of giant or dwarf stars in either cluster may allow us to 
reconcile their copper abundances. As discussed in Sect. [5] a de- 
crease in effective temperature of 300 K in dwarf stars implies a 
decrease by about 0.4 dex. Therefore for the dwarf stars in NGC 
6752 a decrease in effective temperature by about 225 K while 
keeping the temperature of giants constant, would reconcile the 
copper abundances of the to sets of stars. While not implausi- 
ble, this change would certainly cause a mismatch in the abun- 
dances of other elements between giants and dwarfs, most no- 
tably iron, which would become less abundant in dwarfs than in 
giants. While one could argue that atmospheric phenomena such 
as diffusion may alter abundances of dwarf stars, it seems then 
contrived to invoke a different behaviour between copper and 
iron. But let us now turn to the other cluster, NGC 6397. Here 
the situation is reversed, the dwarfs display a lower abundance 
than giants. However here one would need to invoke an increase 
in effective temperatures of the dwarf stars by over 500 K, plac- 
ing them at T e ff around 6700 K and the cluster turn-off at about 
6900 K. While these exceedingly high temperatures may be ap- 
pealing (one would immediately solve the cosmological lithium 
problem!) they appear impossible to reconcile with the colours 
of the cluster and theoretical isochrones. 

Although the precise value of the atmospheric parameters 
assigned to the stars certainly plays a role in the difference in 
copper abundances between dwarfs and giants in the two clus- 
ters, we may dismiss the hypothesis that it may be cancelled by 
a suitable choice of parameters. 

8. Conclusions 

Our study of the Cu i lines of Mult. 1 in the TO stars of the glob- 
ular clusters NGC 6397 and NGC 6752 allow us to draw several 
conclusions: 

1 . the Cu abundance derived from the Cu i lines of Mult. 1 in 
dwarf stars differs from that derived in giants from the Cu i 
lines of Mult. 2 in the same clusters; 

2. the Cu i lines of Mult. 1 in dwarf stars show large 3D correc- 
tions when computed in LTE; 

3. the Cu i lines of Mult. 2 in giant stars show small 3D correc- 
tions when computed in LTE; 

4. the contribution functions of the Cu i lines of Mult. 1 in dwarf 
stars suggest that these cannot be reliably computed under 
the LTE approximation; 

5. for the only star for which we have measured both the lines 
of Mult. 1 and of Mult. 2 we find that the derived abundances 
disagree by 0.5 dex, both using ID and 3D models. 

From the above we conclude that the Cu i lines of Mult. 1 are 
not reliable indicators of copper abundances. A full 3D-NLTE 



treatment of these lines should be used. When treated in ID 
these lines yield abundances which, in the studied cases, differ 
from the abundances in giants by 0.2 to 1 dex, but the differ- 
ence can be in either direction (lower or higher abundance in 
dwarfs). That the dwarf/giant discrepancy is in opposite direc- 
tions for two clusters that differ by only 0.5 dex in metallicities 
also suggests that NLTE alone might not be able to reconcile the 
abundances in the two groups of stars. Whether the combination 
of 3D and NLTE effects may achieve this is an open issue, al- 
though it seems to be doubtful. Clearly a new investigation of Cu 
abundances in giants of NGC 6397 would be highly welcome, to 
place this discussion on firmer grounds. 

Investigation of departures from LTE both in dwarfs and gi- 
ants is needed to place the copper abundances on a firm footing. 
Without knowledge of these departures we recommend that the 
abundances measured in giants from lines of Mult. 2 should be 
preferred for the studies of chemical evolution. The main moti- 
vation for this recommendation is that the 3D effects for the Cu i 
lines used in giants are rather small. 
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Table 5. Atomic data for the Cu i lines used. 
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Table 7. Atmospheric parameters for the gian t stars in NGC 
6752 (lYong et alj |2005|) and NGC 6397 (lGrattonlll982h . 
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Table 6. Copper abundances from the individual lines. 
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Table 8. 3D corrections. 
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